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1 Introduction and Main Results 

Fermion and boson random point processes (fields, or general Cox processus) were stud- 
ied by many authors, in particular since they have a deep connection with the quantum 
statistical mechanics [ini HlHl [TZl El El EEH EES]- See also P [L^ and references 



therein. One of the advantages of the random point field approach to quantum statis- 
tical mechanical models is that it enables probabilistic limit theorems to apply to these 
models. In |ShTa] . typical limit theorems are given for a certain class of random point 
processes which include the particular cases of the fermion as well as boson random 
point processes. In |TI2] . the random point processes, which describe the position dis- 
tribution of constituent particles of boson gases in Bose- Einstein Condensation (BEC) 
are constructed for the first time. 

The purpose of the present paper is to give the limit theorems, such as the Law of 
the Large Numbers (LLN), the Central Limit Theorem (CLT) and the Large Deviation 
Principle (LDP) for the random boson point processes in the regime of the BEC. We 
compare them with the corresponding theorems for the normal phase (i.e. without the 
BEC). In the latter case a detailed study of the limit theorems, which do not use the 
random point processes formalism, is due to [LLS] and [GLMj . In the last reference the 
authors consider even interacting quantum gases, but only in the rarified regime insuring 
the normal phase. These papers motivated the study of the large deviation principle 
in the Bogoliubov-type models [BZj . where BEC plays a key role in description of the 
model thermodynamic behaviour and the spectrum of excitations. 

The study of the boson random point processes in the BEC regime is an interesting 
and delicate mathematical problem [TI2] . see also a recent paper [E]. This last paper 
makes evident that a Cox process in the BEC regime |TI2j is driven by the square norm 
of a shifted Gaussian process. The shift is particle density dependent. In particular, 
this observation makes a contact with the Dynkin isomorphism theorem (known for 
Gaussian processus) as well as a relation between infinite devisability and factorisation 
of the boson Cox process involved in the BEC. 

In the present paper, we study the limit theorems in the BEC regime ( Theorems ll.il 
11.21 ll-3p and discuss in Conclusion the comparison with the analogue of these Theorems 
in the normal phase. 
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Let {G := exp(/?A)}^>o be the (set-adjoint) heat semigroup generated by the Lapla- 
cian acting in L^(M''). For any non- negative bounded measurable function / > with 
a compact support in M.'^ , the operator 

Wf := (G(l - G)-^y/Wl-e~f , 

is a bounded and 

Kf := W}Wf e €i{L'^{M.'^)) , 

i.e., is a trace-class operator on L^(]R'^). If / = 0, then the operator K := Kf=Q is 
bounded with the translation-invariant kerneh 

dp e'^^'-y^P 



K{x,y) 



Below we consider a noninteracting boson random point field Up for the total particle 
density p > pc, which is characterized by the generating functional [TI2] : 

e-W>.,(<ie) = exp(-(.-P.)(Vr^,(l + g/)-'Vr3F7)) ^^^^^ 

!) Det[l + Kf\ 

Here Q{R'^) is the space of all point measures on M , Det stands for the Fredholm 
determinant and (/, ^ = ^ . /(xj), if ^ = Y.j e 



The critical particle density, pc := pdP), for the perfect Bose-gas can be expressed 



as: 



pciP) ■■= K{x,x) 



dp e-^lfl' C{d/2) 



{2ttY 1 - e-/^lpl' (47r/3)'^/2 ' 

The random point field Up defined by by the generating functional (11. ip was intro- 
duced in |TI2] to describe the Bose-Einstein condensation in the non-interacting (per- 
fect) boson gas. For the detailed presentation of these notions, we refer to Ref . |TI2j . 
(See also the next Section O) 

Below in the present paper, we use the following notations: || ■ ||p for Lp{M.'^) norm 
and II ■ II for the bounded operators norm on L'^{M.'^). 

With these notations the main results of the paper can be expressed as follows: 

Theorem 1.1 (Law of Large Numbers) For n — » oo, the limit 
holds m L'^{Q(R'^),Up) . 

Theorem 1.2 (Central Limit Theorem) Let 



v/2(p-p,)||(-/5A)-V2/||2«:(^4-2)/2 

Then the limit 

lim \ e^'^^UpidO = e-*'/' . 
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Theorem 1.3 (Large Deviation Principle) There exists a certain (bona fide) rate 
convex function / : M i-^ [0, +00], such that the limits 

lim sup log ^p(~{f{'/i^)y^)^^^ ^ — inf J(s) for any closed F C M , 



and 



liminf logt/pf — t(/( ■ /k),^) E G) ^ — inf lis) for any open G C M , 

hold. 

In Section 5 we compare these results with those for boson random point processes in 
the normal phase. 



2 Preliminary arguments and general setting 

Let i? be a locally compact Hausdorff space with countable basis, and A be a positive 
Radon measure on R. We suppose that the non- negative (possibly unbounded) self- 
adjoint operator K in L^{R, A) satisfies the following condition K'. 
Condition K' : 

(i) (locally trace class) For every bounded Borel set A G R, K^^'^xa ^ Hilbert-Schmidt 
operator, where xa denotes the multiplication operator corresponding to the indicator 
function of the set A, which we denote by the same symbol. 

(a) The operator G = K{1 + K)~^ has non-negative integral kernel G{x,y), which 
satisfies the conditions: 

G{x, y) > A ® A - a.e. (x, y) e R^ , 

Jj^ G{x, y) \{dy) ^1 A — a.e. x E R . 

The above conditions are arranged in such a way that one can simultaneously deal 
with the random point processes /i^* and ^k,p , see |TI2] and |ShTaj . In particular, the 
operator K has a positive kernel K{x, y), i.e., 

K{x,y)>0 \®\-a.e.{x,y)ER^, 

see |TI2] . The operator K\ := (K^^'^xa)* K^^'^xa is a trace-class operator. For a bounded 
measurable function / with compact support, we define the operator 

Kf := - e-Zi^AV^l - e-/, 

where supp / C A. Note that Kf is independent of the choice of A, which contains 
supp /. 

Let Q{R) be Polish space of all locally finite non-negative integer- valued Borel mea- 
sures on R. Recall that the Borel probability measures on Q{R) (i.e. random point 
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processes on R) /i^^*^ and ^k,p are introduced in |ShTa[ ITI2] for /? > by means of 
generating functionals: 



/ e-'^f'i^ lif'\dO = Det[l + Kf]-' = Det[l + (1 - e-^i^A]"' 
Jq{r) 



(2.1) 



/ e-<^'^Vx,p(ciO = exp{-p(v/l - e-f, (1 + KjyWl - e'f)} 
Jq{r) 

= exp{-p(xA, (1 + (1 - e~f)KAr\l - e^f))} , (2.2) 

for any function / with A D supp /. 

It was shown [IE] that for R = W^ the b oson random point processes corresponding 
to the ideal Bose-gas in the regime of Bose-Einstein condensation (p > pc) is described 
by the convolution Up := fij^^ * fiK,p-pc- 



Theorem 2.1 For any non-negative bounded measurable function f on R with compact 
support supp f G A in a bounded Borel set A G R one has the following equalities: 



(1) 


jQiR) 


= exp[-p(xA, (1 + 


(l-e^O^A)"'(l-e* 


0)], 










'exp[p(\/ e-f — 1 


(1 _ Ve/ - IK^Vef 




-1)]< 


(2) 


[ e<^'«VA',p(c?0 
JQiR) 


= < 


oo 


for WVe^ - 


■ IKAVef - 


1|| < 1 








>. 


for 1 Ve-^ - 


■ IK^Vef - 


1|| ^ 1 


(3) 


JQiR) 


= Det[l + (l-e*0 


Ka]-\ 












'Det[l- Ve/- 


IKaVc^ - 1]-^ < oo 






(4) 


JQiR) 


= < 


oo 


for WVef- 


- IKAVef - 


1|| < 1 



for WVe^^KAVe^^W > 1. 



Proof : Let / ^ 0, i.e., A(supp/) > 0. In |TI2] . pp. 213-214, it was introduced a family 
of symmetric non-negative functions {aA"}n>o defined by the equations: 

exp [-p(x/l-e-/, (l+A7)-Vl-e-/)] = exp [-p(xA, (l + (l-e"^)KA)-i(l-e-^))] 

oo 

= exp [ - p(xA, (1 + i^A)-'xA) + + A^a)-'xa, e-^(i2Ae-^)'(l + i^A)-'xA)] 

1=0 

oo „ 

= / ^A"(a;i,--- ,xOe-^^--^("'=)A®"(dxi---dx„). (2.3) 
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Here, Ra = K\{1 + K/^)~^ satisfies ||-Ra|| < 1 since K/^ is a bounded non-negative oper- 
ator. Using {o"A"}n>0) the random point processes fiK,p was defined as the probability 
measure such that 

f. oo „ n 

/ F{OfiKAdO = y2 ^A"(a;i,--- ,a;„)F(V5,JA^"(da;i,--- ,dx„) (2.4) 

JQ(R) „-n J A" ,_i 



n=0 j=l 



holds for any bounded (or non-negative) measurable functional satisfying F{C,) = F{C,a), 
where ^{A) =^{AnA). 

From this construction, we obtain the first claim (1): 

jQiR) „-n A" 



exp[-p(xA, (1 + Ka)-\k) + p 5^((1 + Aa)-1xa, e*^(/?Ae^O'(l + Ka)-\k)] 



1=0 



= exp[-p(xA, (1 + (1 - e'-f)KAr\l - e^O)]- 
If z G C satisfies |2;|ell-^ll°° ^ 1, then we get the equality: 

oo „ 

J2 / ^A"(xi,-- - ,x„)z"e^"-^(^^)A«'^(rfxi,-- - ,rfx„) 

oo 

= exp[-p(xA,(l + /^A)~'xA)+p5Z^''''((l + /^A)"'xA,e^(^AeO'(l + i^A)~'xA)], (2.5) 

1=0 

Since in the both sides all coefficients the 2;-power series are non-negative, this equality 
( 12. 5p also holds for ^ = 1 in the sense that either the both sides are finite and equal or 
they are both diverge to +oo. When they are finite, we obtain 

/ e<^'«V/<,p(rfO = exp[p{Vef - 1, (1 - Vef - IK We-f - lyWef - 1)] , 

JQiR) 

cf. the proof of Theorem 2.1 in |TI2j . pp. 213-214. Hence, for the second claim (2) it is 
sufficient to show that 

the finite RHS of ([23D ^ \\e^/^RAe^^^\\ < 1 ^ \\^/ef - IK^Ve^ - 1|| < 1 • (2.6) 

Notice that by Proposition 2.3(ii) |TI2] the Condition K'(ii) ensures: RA{x,y) > 0, 
for A ® A-almost all {x,y) G A^. Since R\ is a compact symmetric operator, it follows 
from the variational principle that ||e'^/^i?Ae'^^^|| is the largest eigenvalue of the operator 
C'^^^RaG'^^'^ with eigenfunction (fQ > (A — a.e. on A). Hence we have 

((1 + J^a)"'xa, e^(i?Ae-0'(l + A^a)"'xa) 
= |(y.o, e^/^(l + KAr'xA)\' ||e-^/^/?Ae-^/l' + 0{\\ef/'RAe^/'m 
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for some 6 G (0, 1). Note that | ((/Jq, e'^/2(l + Ka)~^Xa)\ > because (1 + Ka)~^Xa > 
(A — a.e. on A) and ||(1 + K\)~^xa\\ > 0. Thus, we get the first equivalence in fl2.6p . 
For the second equivalence, it is enough to prove that 



WRTe^RTW < 1 ^ \\K]/\e-f - || < 1 

by duality. Let \\R]l^ef R]^"^]] = t] < 1. Then K^^O, f^O and 

1 - Kl/\ef - l)Ki;' = (1 + _ i?f e/i?f )(1 + K^y/' , 

together with R^^"^ R^l"^ > 0, imply 

1 - Kl^\e^ - ^{l + Kj,){l-r])^l-r] . 

Hence K^/^ie^ - l)i^A ^ ^ < 1- On the other hand, if \\K]/'^{ef - l)i^A 1 = ^ < 1, 
then 

1 - R'/'e^R]!' = (1 + Ka)-'^'{1 - K]l\ef - l)J^f )(1 + K^Y'" 

i + ||-f^A|| 

which yields 

< Rfe'Rf « 1 - < 1. 

This finishes the proof of claims (1) and (2) of the Theorem concerning the measure 

The claims (3) and (4) concerning the measure /i^'^*'' can be shown similarly if one 
uses, instead of (12.31) . the representation: 

Det[l + KfY^ = Det[l + /\A]-^Det[l - e'^ Ra]^^ 

CO 

= Det[l + i^A]-' / Per{i?A(x,,Xfc)}i^,-fc^„e-^"-i^(^')A«"(dxi,-- - , 

Til I r>n 

where Det is the Fredholm determinant and Per is the permanent of the corresponding 
matrices |ShTa] . □ 



3 Operators 

Below we deal with the boson random point processes which describe the position dis- 
tribution of the perfect Bose-gas (M"' for d > 2) above the critical particle density 
Pc '■= Pc{f3)i i-e. in the regime of the Bose-Einstein condensation. 

To this end we set R := M.'^ and Rf^ := ^^(1 - C^y^ for K, where := e^^ for 
G. Here /9 > is the inverse temperature and A denotes the d-dimensional self-adjoint 
Laplacian operator in the space L^(M'^) equipped by the Lebesgue measure. Then it can 
be shown that operator satisfies the Condition K', see |TI2j . 
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In the present section, we derive some miscellaneous properties of the operators, 
which we use in the line of reasoning of the next section. First we adopt the following 
definition of the Fourier transformation: 



h{p) := {^h){p) = I e-'P-^'hix) 



dx 



(27r)'^/2 

for h G L^(R'^) and for its extension to L^(R'^). 

Lemma 3.1 For any compact A C M°' the operator (—A)^^^'^xa, {I^^Y^'^Xa is bounded. 
Therefore, 

(-A)- := ((-A)-V2^,)*(-A)-V2xA 
are bounded non-negative self-adjoint operators. 

Proof : These properties can be verified with a help of the Fourier transformation. For 
any g G L'^{W^), we obtain: 



\xAgip)\\ 

12 



^ / / \xrg{p)\'dp 

^ CillxA^lll + ||XA^||2 ^ C2||XA||2lkll2 + II XA || L 11^ II 2 

= (l + c|A|)||^?||i 



A similar argument is valid for the operator K^. □ 



Thus, (-A)-i/2^A is bounded and II (-A)-i/2^a|| ^ y^l + c|A| holds. It gives ||(-A);^^|| ^ 
1 + c|A|. Here, |A| denotes the Lebesgue measure of A. 

'■A- 

Definition 3.2 For k > 0, we define the transformation 

: L\R'^) 3gi-)^ K'^'^g{K-)e L'^iW^). 

Lemma 3.3 The transformation is unitary on L'^iW^) for any k > 0, and it has the 
following properties: 

(1) U^h = h{K ■ ) for the multiplication operator by function h . 

(2) f/.Af/-i = K-^A . 

(3) U^{-/\)~lU-' = k\-A)i\ U^G^U-' = G^'-' . 

(4) U.Kt^U-' = K'J^^ ■ 

Proof : These properties are a straightforward consequence of the relation ^f/^ = U^^,"^ 
on L2(M'^). □ 
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Definition 3.4 For bounded non-negative function f with a compact support and for 
hi > 0, we put 

ft'\x):=±n\e^f^^y-'-l) . 
Lemma 3.5 One has the following estimates: 



X{f>o}{x) 



fix) 



V /(^) 



2^2 



ll/-/l^^l|oo^ 



2^ell/IU/K2_ 



Proof : These estimates are a direct consequence of the elementary inequalities: 



ley - 11 



|e^ — 1 — y| ^ \y\e 



\y\ 



\y\ \y\ 

and \^/z~l\^ for y G M - {0},2 ^ 0. 

Lemma 3.6 For any k > we have the estimates: 

^ (-/?A)X' - ^ (2 



□ 



Proof : Using the Fourier transformation, we get 



p\p\2 ef^\p\'/^' - 1 



\Xa9{p)\'^ dp. 



Then lemma follows from the inequality 



1 1 1 
^ ^ - , for 2/ > , 

y ey -1 2 ' ^ 

and from the estimate ||xa^?||2 = HxAS'lb ^ Hfi'lh • D 

Lemma 3.7 Suppose that supp/ C A. Then for k, ^ one gets the operator-norm 
asymptotics: 



in the space L'^(M.'^). 
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Proof : From Lemma 13.11 13.51 and 13. 6^ we obtain 



^ \\i^/f-^^)i-P^)A'^/f\\ + IIV^(-/9A)x^(a/7- \/^: 



< (llv^lloo + ll\/^l|oo)||(-/3A)x^|| IIa/7- \/^ 



/i^^llLll(-/5A)X^-^-^/^r'll=0(/.-^). □ 



Lemma 3.8 The operator K^^"^ E €i{L'^{W^)), i.e. belongs to the trace-class operators 
on L2(M°'), and 

Tr i^/jK^/^" ^/f] = n^p, [ fix) dx . (3.1) 

Proof : Let { </)„ }„ be a complete ortho-normal system (CONS) functions in L^(R'^) and 
let g{x) = e^^'^XA{x)- Then we have 

}_,\xaUp)\ =1.7^ =(^ = 7^- 

This yields 



n 



Since K^^""^ > 0, it follows that K^^"' G €i(L2(M'^)). Similarly, we obtain the explicit 
value □ 

Lemma 3.9 The operator K^'^ > verifies the following Hilbert- Schmidt norm esti- 
mate from above: 

ll^r^IlL ^ c^nypf'^'^/'il + |log(KV/5)|)|A|(l + |A|) . (3.2) 
Here Cd is a constant depending only on the dimension d > 2. 
Proof: By the Fourier transformation, we obtain 

||r^/3/«^|2 ^ / dq r _dp IXA(p-g)P .ooN 

" "^^ A. (27r)'^/2 y^j, (27r)^/2 (g/^bPA^ - l)(e/3klV-2 - 1) ^'""'^ 
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(27r)'^/2 (27r)^/2 (e/3b+9PA^ - l)(e/3kPA^ - 1) 
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(3.4) 



1°. Case : 2 < d < A. 

From (13 ■4p . we obtain the estimate: 



UT^I^/^^\\2 ^ f dq f dp k,^\xa{p)\ 



^^^^^ '' |Xa(p)P dp f dq 



\p\^~^ {27rYJ^,\e + q\ 



/7 I 2 I xy I 2 



'^dp+ \XK{v)?dp 



'bKi \P\ J\p\>i 



c(||xa||? + ||xa||^) = c(|Ap + |A|) 



Here we changed the variable q = \p\q in the first equahty, and we denote by e a unit 

vector in W^. 

2°. Case : d > A. 

We apply the Cauchy-Schwarz inequality to (13. 3p to get 



|^/3/k2||2 ^ If f \Xa{p - qW dpdq jff \xa{p - qW dpdq 




(e/3blV«^ - l)2(27r)^ ^ J J {emV^' - l)2(27r)^ 

lxX(p)p 



2 dp f dq 



{27rY V {el^liW-' - 1)5 
- -IAl(^) . 

3°. Case : d = 4. 

Let us decompose (13.41) into two parts: 

W^rrns = I l^AXAip)?^ f . f \ dq 



= X1+J2. 

For Ti, |g| ^ 2|p| implies |p + q'| ^ Iq"! — |j5| ^ |q'|/2. Therefore, it follows that 
dp f-^f f dq 



^1 ^ / 777^IXa(p)I 



l^^^ L -^l^kl^2|p| J«;/V;a>|g|>l ^|g|>«/v^ 

dq 



(^g/3|gP/4«2 _ i)(e/3kPA2 _ 1) 
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dp 



\rA{p)\'\9{l-2\p\) 
dq 



Ak^ dq 



Ak^ dq 

?2 1^14 



|xI(p)r^(l-2H)-Cilog 



dp 

V/^v* y|,|>i(el«-|^/^-l)(el«-|^-l) 



2H 



log 



K 



/5 



/3y 



log TTTT + ^ ^2 



bKi/2 



2\p\ 



\(3) 



log 



/k2x4/2, 



^ C4(-j (l + |log(«:V/3)|)(|A| + |Ap) . 



For X2, we obtain: 

T2 ^ 



dp 
dp 



o{j) |A|. 



ixx(p)r 



K'^dq 



\q\<2\p\ /S^b + gPkl^ 

K'^dq 



C3 IIxaII^ 



Thus, we have obtained the desired estimate (13.21) for all cases. 



□ 



4 Limit theorems for BEC 

In this section, we consider the boson random point processus (perfect Bose-gas) in the 
regime condensation, i.e. when 

P > Pci= pciP)) and Up = fi'-^f * ^iKP^ip-p,) , 



where p^P) = K^{x, x) = j^^ ^^i^p _ ^ 



dp 



Proposition 4.1 For a hounded measurable set A C M'^ and non-negative bounded func- 
tion f with supp / C A, one gets the equalities: 



and 



{f\i)^p{di))\p{di) = p [ f{xfdx+Ti[fKifKi]+2{p-p,){f\Kif). 
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Proof : Let us put 

e-^(^) := [ e-<^'«V,(dO , 
then from (I2.ip and (12.20 we get 

= (P - Pc)(XA, (1 + (1 - e--Oi^f)-'(l - e-0) + logDet[l + (1 - e-^i^f] • 
For small e > 0, this yields the expansion: 

W{ef) = ep[ f{x)dx-'\ [ f{xfdx-'\i[fKifKi]-e\p-p,)U.Kif) + 0{€'), 

which implies the proposition. □ 

Corollary 4.2 Under the same conditions as in the Proposition 4-1, one obtains, for 
large n, the following asymptotics : 

and 

= 2«:'^+2(p - p,)(/, (-/3A)XV) + 0(«:^^'^log«:) . 



Proof : Using the unitary operator f/^, we get 

= Tr [/</'^V<'^^] ^ \\f\\lo\\K'J^^\\ls = 0{^'''\ogn) 

and 

= «:^(/, J^f "V) = '^'^'(/, (-/3A)XV) + 0{k''). 

Here we used Lemma [3.91 and Lemma [3.61 Note that supp/( ■ / n) C k,A. Then Proposi- 
tion |4]T] yields the Corollary. □ 

Theorem 4.3 (The law of large number) For n ^ oo and for any bounded func- 
tion f with compact support the limit 

holds m L'^{Q{R'^),Up). 
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Proof : This is a simple consequence of the Corollary 14.21 □ 
Theorem 4.4 (Central Limit Theorem) For k oo the family of random variables 

y/2(p-p,)(/,(-/5A)xV) 
converges in distribution to the standard Gaussian random variable. 
Proof: By Theorem I2.1( l).(3). we obtain 



exp 



[zA«:-('^+2)/2 ((/(■/«:), O-P^:"^ / f{x)dx)] 



exp [ - iXK^'^-^'^/^p [ f{x) dx - W^] , 



where 



W. = {p- p.)(x.A, (1 + (1 - e^^-'^^^^^^/(-/^))K^^j-i(i _ e^A.-(-^)/V(-M))) 

+ logDet[l + (l-e^^-™^(-/^))i^'fj. 

By definition of transformation f/^ and by Lemma 13.6^ the first term can be expanded 
as 



(P - Pe)(t/.X.A, f/.(l + (1 - e^^-^^^^'^^/(-/«))^^j-i(i _ e^A«-(-^)/V(.A))) 
-zAlp-pJ^^-^-^)/^ I fdx + iXK-^'^^y^f^K'J^'f)] +0(1) 



= -zA(p - p.)k^'^^^'^ j fdx + X\p - p,){f, (-/5A)XV) + o{l). 
Here we applied the bound: 

||(i-F)-i-(i + r)IKc||rf 

valid for operators with small enough operator norms with a bound defined by c. 
Similarly, we get also the representation for the second term: 



log Det 



1 + fl -e 



where 



and 



= _,A«:-('^+2)/2Tr[/i^f'^']+^: 
Tr[/Kf'^'] =Pck''J f{x)dx, 



\R\ < Tr 



((1- 
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= 0(A^«:"'^-^)||/||L||i^r'||L = o(l). 
Here we used the bound: 

|logDet[l + F] -TrF| = | logDetsll + F] | = 0(||F||^5) (4.1) 

for the trace-class operators with small operator norms. Recall that Det2[l + Y] := 
e-^''^Det[l + y] = Det[(l-F)e-^] denotes a "regularized" determinant for the Hilbert- 
Schmidt operators Y, see e.g. |ShTa] . 
Thus we get 

= -^Ap«:(^-2)/2 lfdx + X\p - pe) (/, (-/9A)X V) + 0(1) 

and 



exp [^A«-('^+2)/2((/( - Pf^' I , fix) dx)] 

= p-\Hp~Pc){M~PA)-'f)+o{i) 



Then setting A := t/ y 2(p — Pc){f, {—f3'^)\^f), we finally obtain the limit: 

which finishes the proof of the Central Limit Theorem. □ 

Remark 4.5 The above calculations show that the value of the variation that we need 
to normalize the limiting random variable, is contributed from the measure fJ'K^,{p-pc)- 

Before to pass to the Large Deviation Principle, we prove the following lemma. 

Lemma 4.6 Let \\^/f{—pA)]^^^/f\\ < 1. Then — /3A— / is a self-adjoint operator, which 
satisfies the property : Spec {—P A — /) C [0, oo). Moreover, the operator (— /5A — /)^^ 
is bounded and we have: 



(v^, [1 - v^(-/3A)X^/7]-^ //) = / fix) dx + (/, i-pA - f)l'f) 



(4.2) 



Proof : Since the operator —[3 A is self-adjoint, the spectrum Spec (— /3A) C [0, oo) and 
/ is a bounded function, it is obvious that —(3A — / is self-adjoint and (5 — /3A)~^ is 
bounded non- negative operator for arbitrary 5 > 0. Since / ^ and supp / C A, it is 
also obvious that 

^ v7(5 - /5A)-i/7 ^ v7(-/5A)X^/7 . 
Together with the assumption ||\/7(— /?A)X^v7ll < 1; the operator 

oo 

S:=i5- (3A)~' + {6- pA)~'^fY,{^f{5 - f3A)''^)-^{6 - (3A)-' (4.3) 

n=0 
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is a bounded non-negative operator. On the other hand, one can check that 
{S-pA-f)S = I and 8(6 - pA - f) = Ij^^^^^^ , 

which imphes that S = {6 — (3 A — f)~^. Thus, we have —S ^ Spec (— /3A — /), i.e.. 
Spec —/9 A — / C [0, oo). Let {E{\)} be the spectral decomposition of the operator 
{—/3A — /). Then £"(—0) = 0. Moreover, E{0) = holds. Indeed, if one supposes the 
contrary, then there exists a. ip ^ such that 

ifj e E{0)L\M.'^) and {-(3A - /)^ = . 

Thus, we have fip = —^Aip, which implies that 

fip G Ran(-/5A) = Dom(-/3A)-i 

and 

^ = (-/3A)-V^. 

Hence we get //^ = ^/fi-^Ay^fifj = {^{-pA)l^^)^^. This contradicts the 
estimate \\\fJ{—(3A)']^y/f\\ < 1 because y/fip G L^(M'^) belong to the eigenvalue 1 of 
the operator ^/f{—pA)]^^^/f. Therefore, we obtain densely defined non-negative self- 
adjoint operator: 



:-/3A - /) 



-1 



dE{X) 



X 

The boundedness of (— /3A — /)^^ follows from the estimates: 

d{XA(l>, E{X)xa(I)) 

110112=1^0 



-(3A-f)l'\\ = sup / 



X 



^ sup 



sup hm / — — = sup lim(xA0,^X/ 

d>\\2=l "5^0 Jo (> + ^ U\\2=l 



=1 



^ ^ "^l-||v7MA)-v7ll ^ 
To derive equation (14.21) . we exploit the operator (14. 3 p for 5 | 0: 

oo 

(/, [6 -PA- fr'f) = ( v^, J2{^i6 - /3A)-i/7)"^/7) 

n=l 



-(V7, V7) + (V7, (1 - Vf{s - /3A)-i vV)"' V7) 
/dx+(v^,(l- v7(-/9A)x^v7riA/7) , 
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where we used the convergence 



in the operator norm. The latter is a direct consequence of the spectral theorem and 
the dominated convergence theorem. On the other hand, we notice that for 5 J, one 
gets by the monotone convergence theorem the limit: 



(/, (S-pA-f)-'f) = 
d{xAf,E(X)xAf) 



- d{xAf,EiX)xAf) 
S + X 

Therefore, the equality ( 14. 2p is proven. □ 

Theorem 4.7 For any bounded measurable function / ^ with bounded support and 
for any bounded measurable subset A o/M'^ satisfying supp f (Z A we have the following 
limits: 

P{t) := lim ^log / e*'="'<^(-/")'^Vp(c/0 



'pt J^, fix) dx+{p- p,)t^{f, i-PA - tf)l^f) for t e (-00, ||//(-/9A)X^v^iri), 
_oo fortG [||v7(-/5A)X^v7ir\oo) . 

Remark 4.8 (1) If t < \\^{f3A)l^^\\-^, then from Lemma [^T^ we obtain the ex- 
pression for the function P{t) : 

Pit) = Pet f fix) dx+{p- Pc)t{^, [1 - ty^{-(3A)l'y^]-'^) . (4.4) 

(2) By LemmaslEM andl3l\ the operator ^/f {- (3 A)l^ ^/f is a non-negative compact 



operator. Let {v^n} be a CONS of L'^{M?), which consists of the eigenf unctions of this 
operator. We order the corresponding eigenvalues as 

Ai ^ A2 ^ ■ ■ ■ ^ A„ ^ ■ ■ ■ ^ 

Then the Perron- Frobenius theorem yields 

Ai = ||/7(-/5A)X^V^|| and {^,ip,)>0. 
By the remark (1), we obtain: 

P{t)=Pctl f{x)dx+{p-p,)t V '^'f"^'' for t<Ar\ 



17 



which ensures the (essential) smoothness of P: 

P is a C°° function on (—00, A^^^) and lim P{t) = 00 . 

(3) Below we prove that the limits of the function P for the components fiK,p o-nd 
l^^it^^ 0/ the boson random point processes have the following forms: 

PK,,p{t) := lim 4^ log / e*''-'^^^-/^),^^^ (rf^) 



^ (pti^q, [1 - ty7(-/3A)X^v7]-W) for t e (-00, ||v//(-/3A)X^v71|-^) 

\oo fortG [||v7(-/3A)X^y7||-\oo) , 

and 

Pi^fit) := lim log / e--^</(-/^)'«)/.?;Ve) 

'p,t J^,fix)dx fortG(-oo,||v/7(-/5A)X^v7||-i) , 
00 fortG [||v7(-/5A)X^v7l|-\oo) . 



Proof (of Theorem 14.71) : The proof consists of two parts corresponding to t < and 
t > 0. { The case t = is obvious.) 

1°. For t < 0, it is enough to show that 

P(-l) = lim -^log [ e-""'<^( z/p(rfO 



p / fix) dx + {p- p,){f, (-/5A + /)XV). 



To this end notice that from (12.11) and (12. 2p . together with the unitary transformation 
Uk, one obtains the representation: 



P^(U,^1 - e-«-V(-A), U,{l+^l - e-«~V(V-)J^f^Vl - e-«-'/( -A))"! 
^logDet[l + f/.(l - e--^/(-/'^))7rif/-i] 



= -{p- Pc)( V^, (1 + a/^^-M^'^^ V^)-^ V^> 

-^Tr [y^^r V^] - -i,fogDet.[l + ft^^-^Kr] 
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see Definition I3.4[ Then we apply Lemmas 13.51 13.71 to the first term, Lemma 13.81 to the 
second term and Lemma 13.91 with (14.11) to the third term to obtain: 

P(-l) = -(p - pe)(v^, [1 + ^/f{-^3A)l'^r'y^) -pj fix) dx . 

Now it is sufficient to check the identity: 

(/7, [1 + a/7(-/3A)x1 v^]-^a/7) = / fix) dx - if, (-/3A + /)XV) . (4.5) 

Note that the inequahty 

(-/?A + /)X^ ^ (-/?A)^^ 

yields that (— /5A + /)X^ is bounded. Since the operators (e — j3A)~^, (e — /3A + f)~^ 
are bounded and non-negative for any e > 0, we get 

v/7(e - /?A)-i - V^(e - /5A + /)-i = y^ie - [3A)^'y^y^ie -^A + f)^^^' . 
It gives 

v7(6 - /3A)-i^/7 = (1 + //(e - /3A)-i^/7) v^(e - /3A + fT'^f , 
which implies 

1 - (1 + /7(e - /5A)-i v^)-i = v^(6 -PA + f)-'^/f . 
Hence, to verify fl4.5l) . it is enough to prove that 

V^(6-/3A + /)-i/7 ^ V7(-/5A + /)X^v7 weakly, (4.6) 
V^(e-/3A)-i/7 ^ V7(-/5A)x'v7 in norm. (4.7) 

To show (14.61) . let {-E(A)} be the spectral decomposition of — /3A + /. Since 

holds for (f) G L^(]R'^), the dominated convergence theorem yields the limit: 
1(0, v^(-/3A + /)x'/70) - (0, v^(e - /3A + /)x'v^0)| 

/•OO ^1 1 ^ 

= / [j-Y^Jd{^4>,EiX)y^4>)^0. 
To show (14. 7p . we use the Fourier transformation. Put 

II v^(-/3A)x^ - v7(e - /3A)X^ x/7|| 
19 



When d > 4, we obtain that 

n ^ f ^Tfm , ^ f e\Vf^ip)\\ , 
D ^ sup / dp+ sup / — dp 

\m2=iJ\p\<i p \p\ \m2=iJ\p\^i p 
^ ■^(q||/||i + ||/||oo)-o, 

for e ^ 0. When 2 < c? < 4, we get 

D< sup / ^I'^ll- dp 



[ WfWidp 



as e ^ 0. Here we used the bounds ||/70||oo ^ {2Ti)-'^l'^\\f\\]'^\\(j)\\2 and H^lh ^ 
II /II 00^110 II 2, and changed the integral variable p = ^/e/(3p in the latter integral. 
Similarly for d = 4, we obtain the limit: 

D ^ sup / ^"^11^ i, sup / MkCip 

< /• WfWidp e 

^ c||/|K-^iog(i + ^) + -^||/IU-o 

when e ^ 0. 

2°. For t > 0, It is enough to show 

pn ) = h^^' f^""^ dx+{p- p,){f, i-pA - for ||v/7(-/3A)x^/7ll < 1 , 

\oo for ||v7(-/5A)X^v7ll ^1 ■ 

When II v7(~/5^)a^V7|| < 1; then by Lemma \377\ and Lemma 13731 we have 

IIV^^^-'/^r'VTF'll = ||Ve-V(-^)-l At^ Ve-V(7«) _ 1|| < i 

for K large enough, see Definition 13.41 We also use Lemma [3.31 and Theorem 12. II (2), (4) 
to obtain the representation: 



1 



■d-2 



log / e""'<^(-/'')'^>z/p(rfO 
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--^ log Det [1 - f/,Ve--'/(-/-) _ li^f^Ve«''/(V-) - if/^:^] 



+ lTr - log Det, [l - /f ^.^^Af '^^j . 

Applying Lemma 13.71 13.51 to the first term, Lemma 13.81 13.51 to the second term and 
Lemma 13.91 to the third term, we get 

P(l) = (p-p.)(v^, [1 - v7(-/9A)X^v^]-^/7) +Pc / f{x)dx. 

Then the Lemma IT6] proves the case || v7(— /5A)^^v7ll < 1- 

When \\y/f{-pA)l^y/f\\ > 1, we apply and Lemmas [321 EJ] to find that 



||^e/(-/«)/«^ - lA'f^Ve/{-/«)/-^ - 1|| ^ 1 , 
for K large enough. Therefore, we get from Theorem 12. 1( 2) . (4) that 



lim / e<^( '^^^(6/0 = oo . 

gels'*) 



When ||a/7(— /3A)^^v7ll = 1; then applying Lemma [3?71 and transformation U^, we 
find for large k the estimate: 



||v/e/(V«)A^-li^'f^Ve/(V-)A^-l|| = || V/f ^/^-^iff '^V/^^ || ^ 1 - c/^"^ . 

In fact, it is enough to consider the case where the above quantity is smaller than 1. In 
this case Lemmas 13.51 13.71 and 13.11 yield 



+l(v^,{(V^'^-'<'^'v^)"-(v^(-/5A);,^v^)nA/7)l 



^ Mi2^(l + ell/IUA^)eil/ll-/«^(v7, //) 



n + 1 



^ ^ 2 
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This estimate together with Theorem 12.1( 2). give the hmit: 
1 



log 



P- Pc 



d-2 



(v/e/(-/-)A^ - 1, (1 - y^ef - lKf^Ve/(-/«)/«^ - ly' ^/ef - 1> 



= ip- Pc) /^^^ V /^^^) V/^^^) 

n=0 

oo 

> (p-Pc)5^{(A/7,(A/7(-/9A)x^A/7r/7) 



n=0 



n=0 



2c 



oo . 



when K — > oo. Here we apphed in the second equahty, and then the fact that ip is 
the eigenfunction of the operator a//(— /3A)^^-\/7 with the largest eigenvalue 1. Note 
that (vT) V') > 0. In fact, since the integral kernel of this operator is positive on the set 
{/ > 0}, one gets: (p > a.e. on {/ > 0}, c.f. Remark 1131(2). 



The corresponding estimate for n^^f^ is straightforward. 
Recall that the Fenchel-Legendre transformation of the function P has the form: 

I{s) := sup {st - P{t)) 



□ 



By virtue of the Gartner-Ellis theorem (see e.g. |DZj ). we obtained the following large 
deviation principle. 

Theorem 4.9 (Large Deviation Principle) The random variable {f {■/ k) / , sat- 
isfies in the condensation regime p > pc the large deviation principle with a bona fide 
rate function I : 



lim sup —rjr log Up 



and 



lim inf ■ 

K^OO K 



d-2 



^ - inf I(s) 

s&F 



^ - inf I(s) 



forarbitraryclosed F C 



forarbitraryopen G C 



Remark 4.10 Note that contribution of the point processus /i^^*^ to the large deviation 
property is in a sense marginal, since it only shifts the variable s of the rate function I 
see \4-4\) - Taking into account the central limit theorem, we see that the characteristic 
feature of the limit theorems for the ideal boson gas in the presence of the Bose-Einstein 
condensation is reflected by the convolution with a nontrivial component piKf>,p- This 

(det) 

gives Up = p,'j^p * fiKf^,ip-p,)- 
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5 Conclusion 

To compare our results for the case: p > pc (BEC), we would like to mention here the 
corresponding results for the case p < Pc (normal phase without condensation). 

Let us put Kz := 2(7^(1 — zG^)~^ with z G (0, 1), which satisfies p = Kz{x,x) and 
Vp = Pxf^- Then for p < Pc our theorems take the following form, see jLLSl IGLMl IShTaj : 

Theorem 5.1 (The law of large number) For k — oo one has 

l{f{./K),0—P f f{x)dx in L'{Q{R''),Up) . 



Theorem 5.2 (The central limit theorem) For the random variables 

^ ^Kz{x,x) + K^zi^,x)\\fhK^/^' 

one gets the limit: 

hm / e^'^^UpidO = e-*'/' . 



Theorem 5.3 (Large deviation principle) There exists a certain bona fide rate con- 
vex function /' : M 1-^ [0, +oo], such that 

limsup -^logUpf -^(f ( ■ /k),^) G ^ — inf I'(s) foranyclosedF C M 



and 



hm inf 4j log i^p f ^ (/ ( ■ / n) , ^) G ^ - inf /'(s) foranyopenG C M 

hold. 

We may summarize the difference between Theorems 15. 1115. 31 and Theorems ll.mL3l 
as follows. Let 

D. = l(/(-/«:),e), 

be the a random variable corresponding to empirical "density" of particles localised in 
the region of length scale k. 

For the BEC case one gets: 

(i) The random variable converges for k ^ oo to its expectation value m = 
p J^d f{x) dx in mean. 
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(ii) The law of the random variable ~ ^) converges to normal distribution 
as K — > oo. 

(iii) The law of the random variable manifests a large deviation property with pa- 
rameter k'^"^. 

For the normal phase: 

(i) also holds; (ii) holds but for K'^/'^{Di^ — m), instead of K'^^~'^^^'^{Di^ — m); and (iii) holds 
with the order k'^, instead of k'^"^. 

The comparison shows that there are differences in deviation of density fluctuation 
between the BEC and the non-BEC states of ideal boson gases, which reminds the large 
deviation properties for two-phase classical systems, for example lattice spin models, 
see e.g. [P]. The specificity of the BEC is that it is a quantum phase transition with 
particular quantum fluctuations |LePuj . [ZBj . 
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